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The concept of quantum nondemolition (QND) measurement is extended to coherent oscillations in 
an individual two-state system. Such a measurement enables direct observation of intrinsic spectrum 
of these oscillations avoiding the detector-induced dephasing that affects the standard (non-QND) 
measurements. The suggested scheme can be realized in Josephson-junction qubits which combine 
flux and charge dynamics. 



Quantum coherent oscillations in a two-state system 
(qubit) represent the most basic dynamic manifestation 
of quantum coherence between the qubit states. Mo- 
tivated by potential application to quantum computa- 
tion and conceptual interest in macroscopic quan- 
tum phenomena, significant effort is devoted at present 
to attempts to observe and study these oscillations in 
individual "mesoscopic" qubits realized with Josephson- 
junction systems - see, e.g., One of the most direct 
ways of detecting the coherent oscillations in a qubit is 
to monitor them continuously with a weakly-coupled lin- 
ear detector [|] . Spectral density of the detector output 
should exhibit then the spectral line at the oscillation fre- 
quency which contains information about the oscillation 
amplitude and decoherence rate, and has other interest- 
ing features. For example, the absolute intensity of the 
oscillation line demonstrates directly the quantum nature 
of the oscillations. It exceeds by a factor of two intensity 
of the classical harmonic oscillations of the same ampli- 
tude. Quantum mechanics makes larger intensity possi- 
ble by combining harmonic oscillations of the probability 
with discrete jumps of the oscillating variable between 
the two states of the qubit M. 

The spectral line in the detector output, however, does 
not fully represent intrinsic spectral density of the os- 
cillations. In the simplest measurement scheme the de- 
tector measures directly the oscillating coordinate, and 
thus tends to localize it, introducing extra dephasing in 
the dynamics of the oscillations. Such a "backaction de- 
phasing" creates a fundamental limit, equal to 4, for the 
signal-to-noise ratio of the measurement, i.e., the ratio of 
the height of the oscillation line to the output noise of the 
detector. This limitation makes direct measurement of 
the quantum coherent oscillations in an individual qubit 
difficult, and leads to an interesting question whether a 
variant of quantum nondemolition (QND) technique can 
be used to avoid the detector backaction and to overcome 
the limitation on the signal-to-noise ratio. This work sug- 
gests such a QND technique and develops its quantitative 
description. 

QND measurement technique was proposed first for 
detection of weak forces acting on a harmonic oscillator 



in the context of the gravitational- wave antennas JS^U-Ol . 
and was discussed until now in application to measure- 
ments of various realizations of harmonic oscillators - see, 
e.g., [ [Ti"|]f| . Here the concept of a QND measurement is 
extended to the two-state system. In general, the QND 
measurement is realized when a quantum system is cou- 
pled to a measuring detector through an operator (the 
measured observable) that represents at least an approx- 
imate integral of motion. In this case, the backaction 
by the detector which increases fluctuations and uncer- 
tainty in the variables not commuting with the measured 
observable does not couple back into its evolution. Such a 
"decoupling" of backaction makes it possible to measure 
the system continuously without significantly perturbing 
it. 

This discussion implies that specific scheme of the 
QND measurement of a two-state system should depend 
on the main part of the system Hamiltonian. In what fol- 
lows, we consider the case of unbiased two-state system 
that is the most advantageous for the quantum coherent 
oscillations. In the basis of the two states of the oscil- 
lating variable x oc <j z (e.g., charge or flux states in the 
case of charge |||| or flux |^-|^] qubits, respectively) the 
Hamiltonian of the system is then 

H = ~Aa x . (1) 

Here and below cr's denote Pauli matrices, and —A/2 
is the amplitude of tunnel coupling between the basis 
states. The basic idea behind the QND measurement of 
such a system is illustrated in Fig. 1. In the spin- 1/2 
representation of the two-state system ([lj), its dynamics 
can be seen as rotation with frequency A in the y — z 
plane. To perform the QND measurement, the direction 
along which the spin is measured should follow as closely 
as possible the system rotation. This is achieved if the 
measurement direction rotates with frequency Q ~ A. 
(Since the phase of the oscillation of a spin 1/2 can not 
maintain any semiclassical dynamics, there is no question 
of the phase relation between the two rotations.) Thus, 
one can suggest the following Hamiltonian for the QND 
qubit measurement: 
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H = -^Acte - -(cosQta z + smQta y )f + H . (2) 

Here Hq is the Hamiltonian of the detector which is cou- 
pled to the qubit via the force /, and to simplify notations 
all coupling constants are included in the definition of /. 




FIG. 1. Spin representation of the QND measurement of 
the quantum coherent oscillations of a qubit. The oscillations 
are represented as a spin rotation in the z — y plane with fre- 
quency A. QND measurement is realized if the measurement 
frame (dashed lines) rotates with frequency Q ~ A. 

The Hamiltonian (||) of the QND qubit measurement 
is different from that of the non-QND measurement stud- 
ied so far in theory 0,|||| and used in experiments 
j^-fl in the form of the qubit-detector coupling. The cou- 
pling term in ^ should be contrasted with — <r z f /2 in 
the non-QND case when the detector is coupled directly 
to the oscillating variable. Apart from this difference, 
the two situations should be similar, and the detector 
properties in (||) can be taken to be the same as in the 
previously studied non-QND case: the detector is lin- 
ear, with frequency-independent (in the frequency range 
given by the tunnel amplitude A) response coefficient A. 
These assumptions imply that the force / can be viewed 
as random classical 5-correlated variable with the corre- 
lation function 



(f(t + r)f(t)} = 2nS f S(r). 



(3) 



where the average (. . .}o is taken over the detector density 
matrix, and Sf is the constant low-frequency part of the 
spectral density of /, i.e., the detector backaction noise. 
For more detailed discussion of this detector model see 
Jk|. Under the same assumptions, the correlation func- 
tion of the detector response to the oscillations is: 

K(r) = ^Mt)c(t + r)+c(t + r)c(t))}, (4) 

where c = cos flta z + sin Qta y is the operator of qubit- 
detector coupling and the average is now taken both over 
the detector and qubit density matrices. The time depen- 
dence of c's in eq. (0) combines explicit time dependence 
in their definition and time evolution with the Hamilto- 
nian (0). 



To calculate the correlator (Q) we notice that the ex- 
plicit time dependence of the coupling operator c can be 
written as 



c = e ' a z e ' . 



(5) 



This relation follows directly from properties of the Pauli 
matrices and expresses quantitatively the notion of "ro- 
tation of the measurement direction" in Fig. 1. Using 
this relation one can check that the time evolution oper- 
ator S associated with the Hamiltonian (g) has a simple 
form in the rotating measurement frame: 

S(ti,t 2 ) =Texp|-ij^ dt'H(t')^ 



e i£l<r x (t 1 -t 2 )/2 e -iH' {ti-t 2 ) 



(6) 



where H' is an effective Hamiltonian of the system in the 
rotating frame: 



(7) 



Using the fact (demonstrated more explicitly below) 
that the correlator (^) should be independent of the ini- 
tial density matrix p of the qubit, we can take p in the 
simplest form p = 1/2. Equations (|J) and (||) allow us 
then to reduce the correlator (|j) to the following form: 

K(t) = ^Re(a z a z ( T )} , a z {r) = e iH ' T a z e- iH ' T . (8) 

To find the average of the operator cr(r) in (||) over the 
detector backaction noise /, it is convenient to start with 
the Heisenberg equations for a z (r) with the Hamiltonian 
H' . Averaging the resulting equations for the matrix 
elements <xy of a z (r) with the help of the correlator (|^) 
we get: 

<Tii = iS(ai2 - cr 2 i)/2 , a 12 = iS(au - ct 22 )/2 - T<ji 2 , 

(9) 

and (T22 = —tin- Here S = A — J7, and T = nSf is the 
rate of backaction dephasing of the oscillations by the 
detector. 

As the next step, we need to take into account 
environment-induced energy relaxation/dephasing that 
affects the qubit in addition to the detector backaction. 
Assuming that interaction with the environment is weak, 
so that the characteristic relaxation rate is much smaller 
than the tunnel amplitude A, we can simply add the cor- 
responding terms in the equation (^) for <t z (t). It is more 
convenient to do this not directly in the basis of eigen- 
states of the oscillating coordinate used in ([)]), but in 
the basis of energy eigenstates of the qubit. Transform- 
ing eq. @ into the energy basis (in the spin notations, 
the transformation is the 7r/2 rotation around the y axis: 
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02 - ► Ox , a x — > —u z ) and adding the term responsible 
for the environmental relaxation, we obtain the following 
equation for the off-diagonal matrix element of the oper- 
ator s of the oscillating variable (given by o z (t) in the 
original coordinate basis): 



S12 = iSsi2 - T(si2 - s 2 i)/2 - r e si2 . 



(10) 



To find the correlator (g) we need to solve the evolution 
equations for s with the initial conditions S12 = 1 and the 
diagonal elements of s equal to zero. Equations for the 
diagonal elements show then that they remain zero at all 
r, and si 2 (t) (|To| ) determines the correlator completely. 
This means that the real part of the product ss(t) is 
proportional to the unity matrix, and the correlator (pfl is 
independent of initial qubit density matrix. We also note 
that the different form of the terms describing detector- 
and environment-induced dephasing in eq. ( [To| ) is due to 
the QND nature of the detector-qubit coupling, and that 
while the environmental relaxation rate T e was assumed 
to be much smaller than A, both T e and T can be larger 
than the detuning 6. 

Solving eq. ( [l0| ) with the initial condition S12 = 1 and 
substituting the solution in (g) we obtain the correlation 
function of the detector response to the qubit: 



K(r) 
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-(r/2+r e )r 



coshDi 



2D 



-siahDt 



(11) 
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where D = (r 2 /4 — <5 2 ) . The correlation function ( 111 
determines the spectral density (l/7r) J <ItK(t) coscjt of 
the detector output. Including the constant output noise 
S q of the detector, the spectral density is: 



S(u) = S q 



A 2 (r + r e )(r 2 + rr e + (5 2 ) + r e ^ 2 

4tt {uj 2 ~5 2 +T 2 + rr e ) 2 + w 2 (T + 2T e ) 2 ' 

(12) 



If the detuning S is much larger than the dephasing 
rates, the spectral density reduces to the two Lorentzian 
peaks at oj = ±5. In the vicinity of the positive-frequency 
peak, the spectral density can be written as 



S(uj) = S q + 



A 2 



r/2 



8^ ( w - s) 2 + (r e + r/2) 2 



(13) 



This expression shows that for large detuning, the spec- 
tral density is close to the one obtained in the usual, 
non-QND measurement ||[l3| . The only difference is the 
frequency shift of the spectral peak from intrinsic oscil- 
lation frequency A by the frequency f2 of rotation of the 
measurement frame. The main goal of the QND tech- 
nique, avoiding the detector backaction, is not reached 
in this regime. The backaction dephasing broadens the 
oscillation spectral line, and limits the height of the os- 
cillation peak relative to the background set by the de- 
tector output noise. Indeed, eq. ( |l3| ) shows that even for 



T e = 0, the maximum height S max of the oscillation peak 
is S m ax = A 2 /47iT, and the limitation on the response co- 
efficient of the detector from the linear-response theory, 
A < AwSfSq |jl3), shows that the peak is limited as in the 
non-QND measurement: S max /S q < 4. 

The situation, however, changes, if the frequency of ro- 
tation of the measurement frame matches the oscillation 
frequency more accurately, so that S <C T. The oscillation 
line is shifted then to zero frequency, and the lineshape 



is: 



S(co) 



A 2 



r e + 7 



2tt uj 2 + {T e + 1 )- 



, 7 ^ 2 /r. (14) 



Qualitatively, 7 in this equation is the rate of rare jumps 
of spin in the rotating measurement frame between the 
positive and negative measurement direction, and we see 
that the spectral line is now broadened not directly by the 
backaction dephasing but by these rare jumps. The most 
important feature of these jumps is that the rate 7 van- 
ishes together with the detuning <5, and the spectrum (|l4| ) 
of the detector output reproduces then intrinsic linewidth 
of the oscillation unaffected by the backaction dephasing. 
Therefore, "rotating" measurement (|^) with frequency ft 
equal to the tunnel frequency A avoids the detector back- 
action and realizes the QND measurement of the quan- 
tum coherent oscillations in a two-state system. Other 
methods of avoiding the backaction dephasing are based 
on control of the oscillations via the feedback |l7j . These 
methods, however, do not represent real "measurement" 
but rather creation of the oscillations. 

The QND technique is required to measure the intrin- 
sic spectral line of the oscillations in a qubit. Such an 
intrinsic oscillation line is typically found in calculations 
that implicitly assume that the spectral density can be 
measured without disturbance from the detector even in 
quantum mechanics - see, e.g., Jl8[ . As we saw above, in 
the case of mesoscopic qubits this assumption is not ob- 
vious, and non-trivial measurement schemes are required 
to observe intrinsic oscillation spectrum. Without these 
schemes, the detector backaction affects the spectrum at 
least within the linear detection approach appropriate for 
typical mesoscopic detectors. 

In mesoscopic qubits, for which the question of mea- 
surement of an individual qubit is particularly relevant, 
the necessity to measure two non-commuting qubit op- 
erators to implement the QND measurement presents a 
non-trivial requirement since only one (basis-forming) 
observable typically has clear physical realization. A 
system where this requirement can be satisfied is the 
Josephson-junction qubit that combines charge and flux 
dynamics. The simplest version of such a qubit is the 
"charge-controlled SQUID" : two junctions included in a 
superconducting loop with inductance L, with the island 
between the junctions that has small electric capacitance 
C (Tig. 2) JTsj] . External magnetic flux <!> e is applied to 
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the loop, and a gate electrode induces charge q on the 
middle island. 



qubit electrometer 
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FIG. 2. Schematic of the Josephson-j unction qubit struc- 
ture that enables measurements of the two non-commuting 
observables of the qubit, a z and a y , as required in the QND 
Hamiltonian (^). For discussion see text. 

Depending on the parameter values, this system can 
act as a charge or flux qubit. To be specific, we consider 
here the situation when the junction coupling energies 
E\,2 are small in comparison to the island charging en- 
ergy (2e) 2 /2C, and the inductance L is also small. Then, 
dynamics of charge 2en on the middle island dominates 
the system, and for q ~ e is reduced to the two-state 
dynamics with the basis states n = and n = 1. If the 
frequency of the "plasma" oscillations of the flux $ in 
the loop L around the external flux <I> e is much larger 
than the energies of the charge qubit, the oscillations 
effectively decouple from the charge dynamics, and the 
charge part of the system Hamiltonian reduces then to: 



H 



eje - q) 
C 



a z — E + cos ip e <T x + E- sin (p e a y 



(15) 



Here E± = (Ei±E2)/2, Pauli matrices act in the basis of 
the two charge states n = and n = 1, and ip e = e^ e /h. 

The flux part of the system still plays important role, 
since the tunneling of Cooper pairs in the charge qubit 
( |i"5| ) produces current in the loop: 



I = -(/+ sin tp e <r x 



J_ COS (p e CTy) , 



(16) 



where Ij = 2eEj/h are the junction critical currents, 
and I± = (Ii ± h)/^- (The last equation assumes, for 
simplicity, that the junction capacitances are equal.) For 
finite loop inductance L, the current ([l6]) creates small 
variations of the flux $ through the loop that can be 
detected by external system, and by monitoring these 
variations one can measure either a x or a y component 
of the charge qubit (|l5|). The possibility to measure all 
observables of an individual mesoscopic qubit provided 
by the qubit (|l^) can be important for many different 
purposes. 

When $ e = and q = e, the qubit ( |l~5| ) provides the 
necessary elements for the QND measurement discussed 
above. In this case, eq. ( |l5| ) reduces to the unbiased 
Hamiltonian (|l|) with the tunnel amplitude E + , and the 
current / = (7_/2)cr y in the inductance L represents the 



(j y component of the qubit dynamics. The a z component 
can obviously be measured through the charge on the 
middle electrode of the qubit. The final step in realiza- 
tion of the QND Hamiltonian (||) is to convert both mea- 
surement components into one physical form (e.g., charge 
or flux) and apply (with sine- and cosine-modulated cou- 
pling strength) to one detector. One way of achieving 
this is to convert the charge signal (a z ) into the flux 
form by another charge-controlled SQUID operated in 
the regime of Bloch electrometer J2CJ (Fig. 2), when all 
internal frequencies of the electrometer are much larger 
than the frequency E + of the charge signal. The sig- 
nal modulates then the quasistationary critical current 
of the double junction system and, as a result, changes 
the current and the flux in the electrometer loop. Once 
the two signals are in the flux form, the subsequent steps 
can be realized using known Josephson-junction circuits: 
modulated flux transformers and magnetometer. 

To summarize, it is possible to design the QND tech- 
nique for measurement of the quantum coherent oscilla- 
tions in an individual two-state system. The technique 
avoids the detector backaction, and overcomes the lim- 
itation on the signal-to-noise ratio of the measurement 
of the spectral density of the oscillations imposed by the 
backaction dephasing. 
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